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ABSTRACT 


A  method  based  on  the  application  of  the  iodine-131  isotope  was  used 
for  investigation  of  turbulent  diffusion  in  the  sea.  A  special  appara¬ 
tus  constructed  for  the  purpose  in  1956  was  utilized  for  measurements 
conducted  in  the  Black  Sea  on  the  basis  of  which  the  numerical  values 
of  turbulent  diffusion  factors  were  determined.  The  tracers  were  ob¬ 
served  with  the  aid  of  recorders  consisting  of  Geiger  •*  'Her  counters 
which  made  it  possible  to  determine  the  concentratio  £  radioactive 
element  in  a  definite  solid  angle.  The  operation  of  the  apparatus  is 
illustrated  by  graphs  and  the  methods  of  recording;  calculation  and 
verification  of  the  results  are  elucidated  by  mathematical  formulae 
based  on  tests. 
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DETERMINATION  OF  THE  COEFFICIENTS  OF  TURBULENT  DIFFUSION 
IN  THE  SEA  WITH  THE  AID  OF  RADIOACTIVE  INDICATORS 


Radioactive  Isotopes  of  some  elements  find  an  increasingly  wide  appli¬ 
cation  in  various  investigations  conducted  at  the  present  time.  Their 
typical  properties,  depending  upon  the  type  of  investigations  and 
tasks  faced  by  investigators,  are  taken  into  consideration.  On  the 
basis  of  their  prooerties,  the  radioactive  isotopes  are  divided  into 
three  basic  groups.  Depending  upon  the  character  of  investigation 
peculiar  to  the  given  element,  they  are  called  ct,  /6~,  and  7-- radiators . 
The  first  two  groups  of  radiators  emit,  corpuscular  radiation  and,  in 
connection  with  it,  they  appear  to  be  intensely  ionizing  elements. 
Because  of  great  energy  losses  per  unit  travel  distance  (mean  free 
path),  fi~  and  especially  o-particles  are  weakly  penetrating  particles} 
whereas  the  -/-rays  emitted  by  the  various  elements  are  considerably 
more  penetrating  particles  because  of  which  they  make  possible  the 
studies  of  relatively  greater  layers  of  matter.  This  difference  ap¬ 
pears  to  be  one  of  the  main  properties  on  the  basis  of  which  the 
elements  designed  for  the  solution  of  experimental  problems  are  selec¬ 
ted.  A  large  section  of  investigations  is  devoted  to  the  utilization 
of  dependence  of  radiation  absorption  upon  the  properties  of  the  matter 
to  be  investigated.  This  type  of  investigation  makes  it  possible  to 
measure  the  "statistical"  properties  of  matter.  The  method  of  utiliz¬ 
ing  radioactive  elements  as  property  indicators  has  found  a  wide  appli¬ 
cation  when  solving  the  problems  in  which  the  dynamics  of  various 
processes  are  investigated.  Indeed,  the  use  of  indicators  having  the  . 
size  of  molecules  opens  wide  vistas  for  the  study  of  most  complex 
processes.  Because  the  ct  and  p- particles  are  readily  absorbed  by 
matters  and  have  short  travel  routes,  the  greatest  prospects  for  prac¬ 
tical  application  are  associated  with  /-indicators.  The  possibility 
of  observing  an  indicator  at  a  relatively  great  distance  from  the  /6£ 
place  where  the  indicator  is  located,  without  dist orbing  the  medium 
to  be  investigated,  presented  a  basis  for  the  application  of  /-active 
isotopes  to  studies  in  the  turbulent  diffusion  occurring  in  the  sea. 

This  paper  is  deveted  to  the  description  of  the  investigation  method, 
apparatus,  and  the  measurement  data. 


The  Choice  of  Radic active  Indicator 

The  need  for  safety  when  working  with  radioactive  isotopes  lends  a 
specific  character  to  the  work.  In  this  instance  the  investigator 
must  take  into  consideration  not  only  the  requirements  for  optimum 
physical  conditions  during  experiments  but  also  the  safety  regulations. 
In  some  instances  the  former  aspect  must  be  sacrificed  in  favor  of 


TRANS-217 


2/65 


the  latter.  The  main  properties  of  radioactive  isotopes  are  as 
follows:  the  character  of  radiation  (  /J  o,  f  )t  radiation  energy, 
and  the  period  of  half-life  of  the  isotope.  In  our  case,  the  choice 
of  /'-radiator  was  necessary  for  the  purpose  of  observing  the  indicator 
at  considerable  distances  from  the  recorder.  On  the  basis  ci*  the  same 
concepts,  it  appeared  necessary  to  use  the  isotope  having  a  maximum  of 
radiation  energy;  however,  the  selection  of  extremely  great  energy  was 
impeded  by  safety  requirements.  In  addition,  the  curve  characterizing 
the  attenuation  of  V-radiation  in  water  demonstrates  that  with  ener¬ 
gies  exceeding  100  xev  the  dependence  of  the  attenuation  factor  upon 
the  energy  is  considerably  less  pronounced.  This  makes  it  possible  to 
select  isotopes  having  a  moderate  radiation  energy.  In  the  given  case, 
it  is  desirable  to  select  an  indicator  characterized  by  monochromatic 
radiation,  which  may  facilitate  the  interpretation  of  the  resultant 
picture .  Because  the  area  where  the  investigations  of  turbulent  dif¬ 
fusion  are  conducted  cannot  be  isolated,  it  is  necessary  to  select  an 
isotope  having  the  smallest  possible  period  of  half-life.  Experience 
has  demonstrated  that  an  isotope  whose  period  of  half-life  equals  about 
10  days  is  well  suited  for  investigations.  The  radioactive  iodine-131 
isotope  satisfactorily  fulfills  all  these  conditions,  although  it  does 
not  appear  to  be  a  monochromatic  radiator  (Ey  =  0.7,  Mev  *  8.1  days). 

Also  the  question  of  the  quantity  of  the  emitted  isotope  at  each  ex¬ 
periment  is  important.  Preliminary  calculations,  which  have  been  veri¬ 
fied  in  practice,  demonstrate  that  0.1  to  1.0  millicuries  is  an  entirely 
satisfactory  magnitude  for  the  conduction  of  measurements.  At  the 
same  time,  the  concentration  of  isotope  created  by  the  emittance  of 
such  a  quantity  of  radioactive  matter  will  in  a  few  hours  t  scone  safe 
in  open  sea, because  the  process  of  intermixing  diminishes  to  a  certain 
degree  the  magnitude  of  concentration. 


The  Scheme  of  Measurements 

In  1956,  Professor  A.  G.  Kolesnikov  at  the  Department  of  Physics  of 
'  the  Sea  and  Continental  Waters  of  the  Faculty  of  Physics  led  the  con¬ 
struction  of  an  apparatus  making  it  possible  to  determine  the  intensity 
of  turbulent  diffusion  along  the  vertical  in  the  sea.  In  the  summer 
of  1956,  the  apparatus  was  used  for  investigations  conducted  in  the 
Black  Sea  (The  Black  Sea  Branch  of  Marine  Hydrophysics  Institute  of  the 
Academy  of  Sciences  of  USSB). 

The  diffusion  of  the  "cloud"  of  radioactive  liquid  in  the  sea  was  sub¬ 
jected  to  concrete  investigations.  If  a  drop  of  radioactive  liquid  is 
.  introduced  into  the  medium  to  be  investigated  (sea,  in  this  case)  and 
the  diffusion  of  radioactive  matter  is  observed,  it  is  possible  to 
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determine,  on  the  basis  of  these  observations,  the  diffusion  factor  in 
all  directions  by  using  the  method  of  mathematical  analysis. 

The  tracers  were  observed  by  recorders  consisting  of  Geiger-MQller 
counters.  The  structure  of  the  recorders  made  it  possible  to  determine 
the  concentration  of  radioactive  element  in  a  definite  solid  ngle. 

Figure  1  explains  the  test.  The  B  and  C  recorders,  as  well  as  the 
source  of  radioactive  isotopes  A,  were  submerged  into  the  sea  so  that 
they  assumed  a  strictly  definite  position  with  respect  to  one  another. 
In  addition,  the  structure  of  the  instrument  guaranteed  such  a  posi¬ 
tioning  of  the  instrument  in  the  current  that  che  direction  of  the 
current  coincided  with  the  plane  containing  uhe  recorders  and  the 
source  of  the  isotopes.  At  a  definite  time  moment,  the  source  emitted 
a  definite  quantity  of  iodine-131  solution.  After  some  time  (t^)  the 
concentration  of  radioactive  isotope  in  area  G^  began  to  exceed  the 
level  of  the  background  and  the  data  unit  3  recorded  the  change  in  the 
concentration.  As  a  result  of  diffusion,  a  considerable  quantity  of 
radioactive  molecules  entered  the  area  G2  in  time  interval  tp.  Con¬ 
sequently,  the  data  unit  C  also  recorded  the  increase  in  radioactivity 
and  rise  above  the  level  of  background.  Subsequently,  the  radioactive 
solution  continued  to  become  diffused  in  tha  surrounding  medium.  Thus, 
in  the  beginning,  each  of  the  data  units  recorded  the  increase  of 
radioactivity  in  the  corresponding  area,  but  later  each  of  the  units 
recorded  the  fall  of  the  level  of  this  area  to  that  of  the  background. 
At  the  corresponding  time  moments,  t-j  and  t^,  the  data  units  C  and  B 
ceased  correspondingly  to  record  the  excess  of  radioactivity.  The 
graphical  distribution  of  radioactivity  during  the  measurements  is 
presented  in  figure  2. 


m 


The  operation  pattern  of  the  instrument  is  such  that  the  appearance  of 
radioactivity  can  be  recorded  only  when  its  value  in  comparison  with  the 
level  of  the  instrument  background  is  exceeded  two  to  three  times. 
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The  Recording  Apparatus 

Data  Units.  In  order  to  obtain  dependable  results,  the  recorder  must 
have  a  small  "aperture"  (G).  This  requirement  appeared  to  be  one  of 
the  main  conditions  that  had  to  be  observed  when  constructing  the  ap¬ 
paratus*  The  recording  units  (Geiger Sillier  counters)  were  placed  in 
a  steel  container  with  thick  walls  (30  mm).  The  cover  was  made  of 
aluminum  plate  5  mm  thick.  Such  a  structure  made  it  possible  to  /67 
isolate  the  recording  units  from  the  water.  At  the  same  time,  the 
sensitivity  of  the  apparatus  in  the  direction  of  the  aluminum  "window" 
was  considerably  greater  than  the  sensitivity  in  any  other  direction 
due  to  considerable  differences  in  the  absorption  values  of  /-rays  in 
the  walls  made  of  iron  and  aluminum.  A  still  greater  anisotropy  as 
regards  the  sensitivity  would  be  gained  in  this  case  by  the  generally 
known  scheme  of  concurrences.  Regrettably,  the  use  of  Geiger-H tiller 
counters  in  these  conditions  did  not  make  it  possible  to  utilize  the 
scheme*  Indeed,  it  is  known  that  the  sensitivity  of  the  counter  when 
recording  the  /-rays  does  not  exceed  1  %  .  Consequently,  the  sensi¬ 
tivity  of  the  scheme  of  double  concurrence  is  of  the  order  E  *  0.01  %. 
Such  a  small  E-value  requires  the  emittance  of  100  times  greater 
quantity  of  isotopes,  which  is  not  recommended  with  respect  to  safety 
conditions  in  the  work.  Therefore,  the  counters  had  to  be  shielded  by 
a  thick  layer  of  iron  and  they  had  to  be  switched  on  simultaneously  in 
order  to  increase  the  sensitivity  of  counting.  Eight  counters  of  the 
type  MC-7*  were  placed  into  the  body  of  each  recorder.  The  recorders, 
which  contained  in  addition  to  the  counters  one  cathode  repeater  for 
each  channel,  were  connected  by  means  of  a  multicore  cable  with  A 
control  desk  where  the  impulses  obtained  from  counters  were  trans¬ 
formed  in  a  definite  manner.  The  transformed  signal  was  further  trans¬ 
mitted  to  a  loop  oscillograph  J10-U*  which  recorded  variations  in  the 
concentration  of  the  radioactive  element  in  area  G.  The  use  of  the 
cathode  repeater  appeared  to  be  necessary  because  of  a  great  capacity 
of  the  cable  and  a  small  power  of  impulses  coming  from  the  counters . 

The  Unit  of  Impulse  Transformation.  The  impulses  caning  from  the 
cathode  repeaters  of  recorders  were  arranged  according  to  duration  and 
amplitude  with  the  aid  of  cascades  of  locked  multivibrators.  Further, 
the  impulses,  having  been  intensified,  were  transmitted  to  the  inte¬ 
grating  nucleus.  A  loop  oscillograph  had  been  placed  at  the  exit  of 
the  integrating  nucleus.  According  to  the  operation  principle  of  the 
integrating  nucleus,  the  current  passing  through  the  loop  is  in  the 
given  case  proportional  to  the  frequency  of  the  impulses  arriving  at 

*The  Russian  letter  symbols  have  not  been  transliterated  in  order  to 
facilitate  the  identification  of  the  symbols. 


the  entrance.  Because  each  impulse  coming  from  the  data  unit  corre¬ 
sponds  to  the  impulse  formed  in  the  course  of  the  integrating  nucleus, 
the  current  passing  through  the  loop  is  proportional  to  the  frequency 
of  the  impulses  arriving  by  the  cable  from  the  data  unit,  i.  e.  to  the 
concentration  of  jf-active  element  that  is  measured  in  the  G  areas. 
Thus,  the  change  of  current  when  leaving  the  integrating  chain  attests 
to  the  variation  in  the  concentration  of  radioactivity  in  the  sea. 

The  processing  of  such  records  enabled  us  to  determine  the  coefficient 
of  turbulent  diffusion  by  a  method  which  is  discussed  below.  The 
structure  of  the  detector  is  presented  in  figure  3* 


The  Source  of  the  Gamma- Indicator .  When  constructing  the  apparatus, 
which  emitted  at  a  definite  moment  the  needed  amount  of  iodine-131 
isotope,  the  following  conditions  were  taken  into  consideration: 

1)  it  had  to  be  compact,  2)  convenient  for  operation,  and  3)  safe  in 
handling.  The  following  device  was  constructed  with  a  view  to  satis¬ 
fying  these  requirements.  The  source  consisted  of  a  rotating  drum  in 
which  ampoules  containing  the  solution  of  radioactive  isotope  were  set 
on  the  basis  of  radii.  The  axis  of  the  drum  was  firmly  fastened  to 
the  data  units.  The  drum  could  be  set  in  motion  at  any  required  time  /68 
moment.  When  the  drum  was  rotating,  the  ampoules  successively  passed 
near  a  blade  that  would  break  them  open.  Thus  it  was  possible  to 
conduct  16  successive  measurements  without  any  additional  charging  of 
the  apparatus.  The  preliminary  examination  of  the  operation  of  the 
source  was  verified  in  the  laboratory  with  the  use  of  dyes.  The  result 
was  completely  satisfactory,  which  made  it  possible  to  apply  the 
method  at  the  time  when  the  radioactive  isotope  was  emitted. 


The  Theory  of  the  Method 

When  processing  the  data,  it  is  necessary  to  obtain  a  mathematical 
definition  which  in  a  satisfactory  measure  reflects  correctly  the  true 
physical  pattern  of  the  isotope  diffusion  in  the  sea.  The  mathematical 
physics  includes  a  problem  on  the  diffusion  of  clouds  formed  as  a 
result  of  explosion  of  charges.  The  difference  in  the  given  case  is 
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characterized  only  by  two  peculiarities.  The  cloud  itself  is  a  source 
of  "luminescence, "  whereas  m  the  problem  that  is  discussed  the  cloud 
absorbs  the  light  coming  from  infinity.  The  second  peculiarity  is 
that,  when  solving  the  problem,  it  is  necessary  to  account  for  the 
absorption  of  light  in  the  medium  between  the  cloud  anc  the  observer  or 
the  instrument  that  records  the  intensity  if  its  luminescence  in  a 
given  direction.  Evidently,  these  peculiarities  are  not  essential  and 
therefore  the  conclusions  can  be  drawn  in  a  quite  analogous  way. 


This  process  can  be  presented  in  the  form  of  the  diffusion  of  a  cloud, 
ignoring  the  gravitational  force  of  the  earth-  and  inequalities  in  the 
initial  scattering.  In  addition,  it  is  assumed  that  the  dependence  of 
the  intensity  of  the  process  upon  the  decomposition  during  the  measure¬ 
ment  can  be  neglected.  In  these  conditions,  the  concentration  of 
radioactivity  in  a  cloud  is  expressed  by  the  following  formula 


x'+y***' 

4DI 


where  Q  denotes  the  quantity  of  radioactive  isotope  emitted  during  a 
single  test,  D  is  the  coefficient  of  turbulent  diffusion.  Thus,  the 
problem  of  diffusion  from  a  momentary  source  of  current  with  power  Q 
in  unlimited  space  arises. 


Let  us  direct  axis  z  along  the  axis  of  vision  and  assume  t>-at  the 
observer  is  found  in  infinity.  In  this  caso  the  cloud  v’-'.xl  be  as  if 
projected  on  plane  x,y.  Let  us  appraise  the  visibility  of  various 
areas  of  the  cloud.  It  is  known  that  the  intensity  of  radiation  at  a 
distance  z  f rom_  the  source  in  an  absorbing  medium  is  expressed  by 
formula  D  *  D0e  ^  ,  where  D0  is  the  intensity  of  radiation  at  point 
z  *  0,  and  /*  is  the  attenuation  factor  of  the  current  of  gamma-rays 
in  water.  If  the  distribution  of  concentration  of  radioactive  mole¬ 
cules  is  given  by  formula  U(x,y,z,t),  the  total  intensity  along  the 
axis  of  vision  is  determined  try 


,  y, t)e~'+*dz=  —~ 


Q_ 

(2  Y  T~Dt)* 


t1  +CO 

Wt  ['  . 


The  integral 


e~ ' 


02 


t  +  2Dfy. 

2  VDi 


is  reduced  to  the  Gauss  integral  by 


introducing 
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Finally  we  have 


A  = 


*Dt  : 


(2  y  xDt)' 


>l  (2  yrut)* 


.  e  4DI  gDtr! 


*2 yot  jY°’ 


'v'  dv  • 


Q  ~  4DI 


4  zDt 


e 


Thus 


4«0< 


e-  •>’  e 


40/ 


where  p  -  *2  +  y2. 

In  the  calculations  we  aid  not  account  for  the  current  speed  because  in 
tne  area  where  the  investigations  were  contacted,  the  phenomenon  could 
be  ignored. 

The  "visibility*'  of  the  cloud  of  radioactive  particles  is  determined  by 
the  value  of  the  threshold  of  sensitivity  of  the  recording  instrument, 
and  thus  the  resultant  formula  is  reduced  to  a  form  making  it  possible 
to  apply  the  'orraila  to  the  processing  of  the  data  that  are  obtained. 
Thus,  taking  into  consideration  the  threshold  of  sensitivity  of  the 
instrument,  the  observed  cloud  has  finite  sizes  at  each  time  moment  and 
its  radius  is  determined  by  the  given  formula.  Finally,  the  equation 
relative  to  the  coefficient  of  turbulent  diffusion  assumes  the 
following  form: 

40®<y  -  4Dt  (in  =  p\ 


Q 

where  a  —  — — . 

4i» 

An  analysis  of  the  Gauss  integral  demonstrates  that  if  the  error  is 
limited  to  10£  it  can  be  assumed  that  the  observer  is  in  infinity  at 
v  >1.17.  If  a  greater  accuracy  is  needed  and  if  the  magnitude  of 
argument  is  smaller  than  the  given  value,  the  method  of  solution  re¬ 
mains  the  same;  however,  in  lieu  of  integral  VT~  ,  it  is  necessary  to 
use  tables  in  order  to  arrive  at  the  correction  formula  for  the 
solution  of  this  problem. 

As  is  seen,  the  final  formula  includes  a  magnitude  that  characterizes 
the  ratio  of  sensitivity  of  the  instrument  to  the  quantity  of  the 
isotope  that  is  emitted.  In  order  to  eliminate  this  magnitude  from  the 
final  solution,  corresponding  measurements  have  to  be  carried  out. 
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Indeed^  let  us  .solve  the  problem  with  the  presence  of  two  radii  along 
which  the  recording  instruments  extend.  Suppose  the  time  moments  are 
tj  and  t2  between  the  emittance  of  isotope  and  the  beginning  of  re¬ 
cording  of  radioactivity  tliat  exceeds  the  threshold  for  each  of  the 
distances  ^  and  respectively.  In  this  case  we  have  the  following 
system  of  equations: 

4D*/y  -  ADtl  In  ■—  =  £«* 

4D*/|u*  _  ADt.  In  — *  =  o?. 

*  a  1 


The  second  equation  is  transfonned  and  ve  have 

AD'&f  -  4 Dtt  In  —  —  ADtt  In  ±  -  $ 
*  *  h 

by  changing 

In  —  —  In  —  —  =  In  —  -f  In  — . 
a  a  tt  a  t, 

Solving  the  system  jointly,  we  find  that 


Hence 


-4D*2ln-^  +  (-^P?-pi)  =  0. 

+  \/^  16^  In*  ~  —  -  Mt)  Pi  —  Pjj 


D  _  B  +  V&-AC 
A  * 


A  =  4|i*  (/|  —  /,/,). 


/70 


The  Results  and  Conclusions 

On  the  basis  of  the  der3_ved  formula  we  processed  the  records  obtained 
as  a  result  of  experiment.  Because  a  series  ol  measurements  was 
carried  out  with  a  view  to  verifying  th9  correctness  of  results  during 
small  time  intervals,  it  was  possible  to  compare  the  results  of  the 
entire  series.  This  enabled  uo  to  eliminate  to  a  degree  the  effect  of 
the  introduction  of  the  isotope  into  the  given  medium.  The  results 


TRANS-217 


9/10 


differed  very  little  throughout  the  series  of  measurements.  The  coef¬ 
ficient  of  turbulent  diffusion  fluctuated  from  one  measurement  to  the 
other,  the  values  lying  between  5  cm^/sec  and  7  cxfi/  sec. 

The  work  demonstrated  that  it  is  possible  to  measure  the  coefficient 
of  turbulent  diffusion  with  the  aid  of  radioactive  isotopes.  In  con¬ 
junction  with  other  methods,  this  method  makes  it  possible  to  visualize 
better  the  pattern  occurring  in  the  sea  at  turbulent  intermixing.  At 
the  same  time  one  asserts  that  tne  apparatus  and  means  needed  for  safe 
work  are  sufficiently  simple.  All  this  makes  it  possible  to  utilize 
the  apparatus  in  expeditionary  conditions.  At  the  same  time,  the  use 
of  scintillation  counters  in  recorders  ir.  conjunction  with  the  schemes 
of  concurrence  will  enable  us  to  make  the  instrument  lighter  and  to 
increase  its  accuracy  and,  consequently,  the  accuracy  of  data.  In 
conclusion,  the  author  expresses  his  gratitude  to  Professor  A.  G. 
Kolesnikov  for  the  assistance  rendered  in  the  operation  of  the  instru¬ 
ment  and  advice  that  greatly  aided  the  processing  of  resultant  data. 
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